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Growth dynamic of real networks because of emerging complexities is an open and interesting 
question. Indeed it is not realistic to ignore history impact on the current events. The mystery 
behind that complexity could be in the role of history in some how. To regard this point, the average 
effect of history has been included by a kernel function in differential equation of Barabasi Albert 
(BA) model . This approach leads to a fractional order BA differential equation as a generalization 
of BA model. As opposed to unlimited growth for degree of nodes, our results show that over 
time the memory impact will cause a decay for degrees. This gives a higher chance to younger 
members for turning to a hub. In fact in a real network, there are two competitive processes. On 
one hand, based on preferential attachment mechanism nodes with higher degree are more likely 
to absorb links. On the other hand, node history through aging process prevents new connections. 

Our hndings from simulating a network grown by considering these effects also from studying a 
real network of collaboration between Hollywood movie actors conforms the results and significant 
effects of history and time on dynamic. 


I. INTRODUCTION 


One of the widely studied model to describe network 
dynamics, Barabdsi- Albert (BAT is based on prefer¬ 
ential attachment mechanism According to this 

model, over time new nodes join system and link to the 
earlier nodes, with probability proportional to the de¬ 
gree of preexisting nodes. Nodes with higher degree have 
greater chance to attract more new connections. This 
leads to the emergence of hubs with very large degree. 
Therefore, BA model characterizes one of the important 
features in real networks, namely, scale invariance 

BA model has been very successful in describing many 
properties of real networks [l|, however, its unlimited 
growth prediction for degree of nodes is not always in line 
with reality ©[HI- In fact, in real networks other fac¬ 
tors work against a node’s growth which BA model fails 
to account for [T2l - [l3| . For instance, in realistic networks 
like citation and scientific collaboration network [if 


world-wide web network 17| and network of movie actors 
collaboration [T^ . there are some important phenomena 
such as aging, screening and censorship to be included 


In network of movie actors for example, overtime, hubs 
are replaced by new ones. Superstars get promoted for 
some time period as they get to their height of popularity 
and eventually loose their attractiveness in eyes of me¬ 
dia and decay gradually. New generation of stars will re¬ 
place the old generation and repeat the same cycle again. 
An other example is political network [2^, [2^ , influence 
of people in politics grow and eventually decreases over 
time, political hubs are replaced by other people. No one 
experience unbound growth in politics. 

As opposed to many realistic situations where aging 
process causes the hubs to be replaced by newcomers, in 
BA model an emergent hub will always remain power¬ 
ful and prevent newcomers from becoming strong hub. 
In other word, we need to consider the effect of losing 


power for old hubs and emergence of new hubs. Number 
of studies have been done in this area and some strate¬ 
gies have been proposed Elavil. As a common issue 
with most of these approaches, the outcome has been 
requested manually from the model, i.e, it is imposed 
on the network by adding some new terms to BA equa¬ 
tion. The objection against these approaches is that, new 
terms that they added are not capable to demonstrate the 
main issue; even more, it is not always possible to find a 
closed form solution. 

According to BA model nodes can connect to each 
other without any restriction, hence their degree can 
grow boundlessly. However, in many real systems, each 
node has limited capacity for joining other nodes. For 
example number of active friends a person can have is 
limited as a person would denote part of his/her time 
to each of those friends. Since we have limited amount 
of time per day/weak, one can not simply have infinite 
friends. 

A scientist or actor has just enough time to collaborate 
with small number of partners. Therefore, having some 
connections makes the space and opportunity scarce for 
more connections. It’s current connections do not allow 
the node to join other nodes easily. This is similar to the 
phenomenon in fermi systems known as screening effect, 
reduction in effective electronic potential because of the 
cloud of electrons around one electron [s^. To extend 
the preferential attachment mechanism, the model has 
to also pay attention to these kind of effects. 

A good candidate to include history effects could be 
fractional calculus. Fractional calculus, the generalized 
form of ordinary differentiation and integration to non¬ 
integer order has unique features such as nonlo¬ 

cality and memory, making it hi ghly ap plicable in many 
fields of science and engineering (Sbi - I^ . With respect to 
the presence of kernel for history in fractional operators, 
the results from this model carries the system memory. 
In other word, memory of the system plays a substan- 
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FIG. 1: Schematic diagram represents dynamical growth of nodes in a network, developed under preferential attachment 
mechanism by considering limitations imposed by aging process and screening effect on growth process in (fTalfT^ and 
comparison with standard preferential attachment without considering these restrictions dHIIfl). Red to bine color shades of 
nodes represent their attractiveness for new links. Highest attractiveness is represented by red. The radial distance from 
central node shows arrival time for a node. Over the time a new node joins the network and connects to m existed nodes 
according to their degree. However, old members because of their age have lower chance of being selected. As show, 

in the absence of these effects hubs are dominated in absorbing new links. 


tial role in determining the future path. The proposed 
approach in this work to generalize BA model is based 
on fractional calculus. In this approach, the equation 
governing the dynamic of network growth is a non in¬ 
teger order differential equation. Hence, operators from 
fractional calculus appear. 

By working in fractional realm, the time distances are 
no longer identical. In fact in the growth process, some 
nodes can face long time delays before jumping to the 
next time step. In other words, they become freezen in 
time. At the same time some other nodes experience very 
short waiting period and pass more time steps relative to 
the frozen ones. 

A simulated network based on fractional order BA 
differential equation approach has also been developed 
in which an old node looses effectiveness of its degree. 
Therefore its attractiveness to absorb new links will de¬ 
crease over time. Figure m schematically represents dy¬ 
namical growth of a network evolved based on preferen¬ 
tial attachment mechanism. Moreover, aging and screen¬ 
ing effects have been considered. As a node gets older, 
its probability to attain more links from new members 


decreases. Existing hubs will decay eventually and new 
hubs can emerge. 

The rest of this paper is organized as the following. 
In section m fractional version of BA model by use of 
fractional calculus operators is introduced. Section m 
devotes to solving the fractional equation and it’s results. 
Results of a network simulation by considering aging ef¬ 
fects on it’s dynamic has been presented and discussed 
in section m- In section m the real network of Hol¬ 
lywood movie actor collaboration which shows bounded 
growth for collaboration rate is investigated. Section m 
concludes the results. 


II. FRACTIONAL APPROACH FOR 
PREFERENTIAL ATTACHMENT ALGORITHM 

In the BA model for network growth, the governing 
equation system, 

dki{t) _ mki{t) 

d-t ~ EpTb) 

hiUo) = TO, 


( 1 ) 
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is an integer order equation for each node i, does not 
include the past history of the considered node. Only 
node’s degree in previous step affects its present degree. 
In order to consider the history on the growth process, 
a kernel on time, nit — t'), will be imposed on the right 
hand side of the Eq. (HD- This kernel accounts for the 
memory [iol - li^ . 


dki{t) 

dt 



dt' K,{t — t') 


mki(t') 


( 2 ) 


In case of memory less systems, kernel should be a Dirac 
delta function, — which results in standard differen¬ 
tial equation in BA model for network growth. Here, we 
choose a power law functional form for kernel on time, 
capable of involving the past. This kernel, guarantees 
existence of scalin g fe ature as it is an intrinsic nature of 
most phenomena |44l - l4q . A distant past event should 
have far less effect on present, compared to near past 
events, only exceptions are large influential events which 
even over shadow recent ordinary events. Substituting 
power law kernel in Q, it becomes: 


dki{t) 

dt 


rli'd 




From fractional calculus it is apparent that the right hand 
side of this equation is a fractional integral of order (a — 
I), 01^ the interval [to,t] [s^. Therefore, it can 

be shown as. 


dki(t'j _ 1 ) 


mki(t) 


(4) 


Now applying a fractional Caputo derivative of order 
(a — 1) on both side of the above equation, we can 
write it in the form of a differential equation. 


C 


D? Mt)] 


mki(t) 


(5) 


III. NUMERICAL RESULTS OF FRACTIONAL 
ORDER GROWTH EQUATION 


In attempting to deploy fractional calculus in BA 
model of growing networks, we start by dynamic equation 
in form of. 


ki{tio) = m. 


kijt) 

E'fcA*)’ 


( 8 ) 


Where 0 < a < 1. For a = 1 the above equation be¬ 
comes the well-known dynamic equation in BA model. 
We have a system of fractional order differential equa¬ 
tions coupled by the summation in denominator. From 
here the problem becomes an initial value problem for 
FDE, 


toOtvit) = 
y{tO) = yo, 


which can be solved numerically by the predictor- 
corrector algorithm Hence, we can reformulated 

it to equivalent Volterra integral equation in the form, 

= + (t-s)^°‘-^^f{s,y{s))ds. (10) 

r(a) Jto 

To deal with this integral, we use product rectangle 
method, which divides the domain into n fragments, 
tj = to+jh, with equal space h. The right hand side func¬ 
tion, in terms of numerical approximation yj for y{tj), is 
denoted by f{tj,yj). Finally, we have the discrete form 
as follows, 

yn = yo + (ii) 

with bn coefficients as. 


b 


n 


{n + 1)“ — (n)“ 

r(a -k 1) 


( 12 ) 


As a result, the discrete form of fractional order differ¬ 
ential equation ([5]) becomes. 


For a continuous function / on the [a, b] inter val, left 
Caputo derivative of order a is defined as follow p^ : 

ID? [fit)] = 

where n is the smallest integer greater than or equal to 
a, n = [a] + l. Caputo derivative has the advantage that 
in solving fractional differential equations (FDE), uses 
integer order boundary or initial conditions. 

In the last step to obtain Eq. we applied the fact 
that Caputo fractional derivative and fractional integral 
are inverse operators for a > 0, 

lD?[t„D?<^]fit) = fit). (7) 

In this way, the governing equation is a fractional order 
differential equation guarantees the presence of memory. 


kn = ko + (13) 

l^j kj[t) 

Here, bnS are time dependent coefficients which in¬ 
dicate the aging effect. This factor shows contribution 
of the previous degree of ith node on its present value. 
With increasing the lifetime (increasing n), becomes 
smaller. In other words, over the time old links of a node 
will loose their effect on its growth process based on 6„’s 
coefficients. Therefore, the effective degree of a node (Eq. 
(f01) l will decrease. For a = 1 the bn converges to unity 
and we get back to the standard BA model in which old 
degrees all have the same weight. By numerically solving 
equation system (fT^ . we found interesting results con¬ 
firming the effect of history on network evolution. The 
results in Fig. (I2ap for various amounts of a, show time 
dependence of ki{t) (effective degree has been left for ith 

















4 




t t 

(a) (b) 

FIG. 2: Figure (I2all is numerical solution of Eq. (US and (Hg is mean degree of nodes from simnlating a network that it’s 
growth dynamic is affected by aging process. It is clear from results that k(t) behaves differently from unlimited growth 
predicted by BA model. Dne to aging process k(t) reaches a peak and declines gradually afterwards. Simulation has been 
done for 5000 time step with initial mo = 10 nodes and every new node connects to m = 10 earlier nodes. 


node at time t). It can be seen that for all values of 
a < 1, effective degree of node i increases at first. It 
reaches the maximum value, then starts to decay. 

This behavior could be as a result of competition be¬ 
tween two mechanisms governing the network dynamic. 
On one hand, nodes absorb new links according to pref¬ 
erential attachment algorithm, on the other hands, aging 
process will reduce the probability of being selected by 
new members by reducing its effective degree at that mo¬ 
ment. The older the member, the smaller is its portion 
from past history. However, on the same time period, 
node with higher degree would receive less impact from 
aging than node with smaller degree. 

Many real world networks exhibit the above mentioned 
behavior. Here as a well known example which will be 
studied in particular in section (fVT) , network of movie ac¬ 
tors collaboration could be mentioned. In this network, 
we can see that superstars do collaborate with many ac¬ 
tors as they become popular. However, because of screen¬ 
ing effect one can just have the ability to cooperate with 
finite number of actors. There is not enough room for 
collaboration with every one. The rate of collaboration 
will reach a peak and then decline gradually over time as 
a result of aging process. 

The growth and decaying rates depend on the order of 
fractional equation. Smaller a reaches the peak quickly 
but larger a takes more time. The smaller the exponent 
is, the more aging can reduce growth rate. Getting closer 
to exponent 1 aging effect is reduced and nodes have 
longer growth time. For a = 1 which results in preferen¬ 
tial attachment without aging, it increases boundlessly. 
Exponential increase in time distance that effective de¬ 
gree takes to reach its maximum for different a has been 
shown in Fig. 

Consequently, the fractional order BA equation for 
growing networks can demonstrate the existence of his¬ 


tory. Because fractional order derivatives by applying a 
kernel over time involve the effect of elapsed time, it indi¬ 
cates the fact that every member has an end and will be 
isolated. In many real networks, members do not remain 
in the system for ever. As time passes, they will be set 
aside from the community. The results from fractional 
equation (fOl) remarks this fact. 


IV. SIMULATING AN AGED NETWORK 


Along side solving fractional order dynamic equation 
governing the network, we have simulated a network with 


10 '' 
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FIG. 3: The time to reach maximum for different 
0 < a < 1 values with the fitted curve. This time decreases 
exponentially with increasing the exponent. By considering 
even a very small effect of memory on evolution, which 
means any deviation from 0 = 1, there would be a peak 
after that degree of nodes will experience a decay. However, 
for closer values to 1 this time tends to infinity. 
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(d) a = 0.3 


(e) a = 0.7 


(i)a = l 


FIG. 4: Mall shows probability distribution function for different orders of fractional derivative. While for a = 1, BA model, 
the expected power law behaviour is found, for fractional orders less than unity, there is deviation from power law. This 
deviation gets severe for smaller a. It could be caused by two competitive processes, namely preferential attachment 
mechanism and screening effect. (|4bl) is closeness centrality averaged over all nodes; comparison between three exponents 
a = .3, 0.7,1. Horizontal axis is nodes according to their arrival time. For higher values of a, this measure has higher 
amounts. Elders have shortest path to others. (14cII indicates clustering coefficient averaged over all nodes for three exponents 
a = .3,0.7,1. Horizontal axis shows nodes according to their arrival time. As it is clear, it has certain behavior. Older nodes 
have more communities around. (I4dl4fll show degree-degree correlation averaged over all nodes for three exponents 
a = .3, 0.7,1. Tendency to bind with similar members have been increased by including history effects. The fact is that the 
inclination to powerfuls is reduced. Powerful sites would not be in power for all the time. There is a competition by new 
generations. Hence, it is not the priority anymore to connect to the elders. 


aging process. In this simulation we start by a fully con¬ 
nected network with m-|-l members, hence each node has 
initial degree of m. In each time step i, node i is intro¬ 
duced to the network, which is linked to m previous nodes 
proportional to their degree. To include aging process in 
the growth dynamic, the effectiveness of old links for a 
node is lessened. A link between two nodes, is less ef¬ 
fective for old node, however it preserves its influence for 
younger one; similar to a connection between two nodes 
in a directed weighted network. As a consequence, the 
probability of being selected by newly arrived members 
will be reduced for members with more old links. This re¬ 
duced effectiveness of older links in growth dynamic has 
been applied by using coefficients which are propor¬ 
tional to link’s life time as weight factors in calculating 
probability of receiving new links. 

Figure (l2bl) displays time dependence of effective de¬ 
gree on time steps resulted from simulation. The ob¬ 
served behavior is in good agrement with numerical re¬ 
sults, (Fig. ((^a1) l. obtained from fractional order differen¬ 


tial equation (1131) . Just like numerical solution, initially, 
effective degree per time of a node increases, after reach¬ 
ing its maximum value it decays gradually. 

In figure (|3]) , we have shown some statistical proprieties 
of the network simulation with aging concept along side 
BA model. Small a represents severe aging effects, there¬ 
fore behavior of quantities with smaller a shows greater 
deviation from BA model. 

Degree distribution reveals major characteristics of a 
model. For preferential attachment it follows a power law 
form [I| , however in presence of aging it shows deviation 
from power law. This deviation from power law could 
be as a result of competition between aging process and 
preferential attachment mechanism. According to pref¬ 
erential attachment nodes with high degree absorb most 
of new links, however, getting older over the time reduces 
probability of high degree nodes to be selected [ll| , since 
aging reduces the effective degree of nodes. This reduces 
growth rate of older nodes which gives nodes with smaller 
degree a higher chance to receive new links. 
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FIG. 5: A sample of total collaboration per year for Oscar winner actors. Collaboration rate is number of staring 
actors/actress a actor plays with in a movie each year. It is clear that each actor would experience a peak and gradual decay 
in their total collaboration as they age professionally. 


The average shortest path to all nodes in a network 
known as closeness centrality is a criteria to measure 
power of centrality for nodes in network. This quan¬ 
tity has been studied for different orders of fractional BA 
equation. It can be seen from Fig. dlEI that closeness 
centrality for nodes has raised for aged system compared 
to BA model. This figure shows that older nodes in BA 
model have more central position in network respect to 
olders in aged system. Moreover, in an aged system, older 
members are more close to others than younger nodes. 


Aging process also has impacts on clustering coefficient 
of nodes. Fig. (H^ . It is an indicator that displays ten¬ 
dency of a node’s neighbors to connect to each other [T^. 
Senior members are associated with communities which 
have high degree of connectivity as opposed to recently 
added members. Although aged nodes lose their impor¬ 
tance, they still preserve their strategic position among 
others. However, in BA model older nodes have smaller 
clustering coefficient than olders in system by including 
aging effect. 


One significant consequence of considering aging is 
that new members can have the opportunity to grow 
and become a hub, as opposed to BA model where only 
early members have a chance of fast growth. Hence ten¬ 
dency of nodes to connect to their similar nodes, what is 
called assortativity will increase. This can be seen 
in degree-degree correlation palette in Figs. (I4dll4fl) . De¬ 
spite BA model that considerable connections are with 
hubs or older nodes, in aged networks remarkable corre¬ 
lation can be found between similar nodes in degree. 


V. NETWORK OF OSCAR WINNERS 
COLLABORATION 

In order to verify the above results for a real network, 
collaboration network of Hollywood movie actors have 
been studied. We extracted list of actors from the In¬ 
ternet Movie Database available at www.imdb.com. The 
table of all movies and star actors of those movies have 
been collected. TV series and TV movies are not in¬ 
cluded in our data. Nodes in the intended network are 
actors and they have a common link if they have acted 
in the same movie. As a case study in this network, to¬ 
tal number of collaboration of Oscar winners with other 
actors per year has been derived as shown in figure (O 
for some of them. It is the total number of staring ac¬ 
tors that an actor/actress played with in each year. It 
can be seen that almost all the actors presented in the 
above network have the same pattern in their career. At 
start of acting career, collaboration rate is increased till 
it reaches a maximum. After that, collaboration rate 
will decrease gradually. Despite some exceptions, it is a 
general behavior that collaboration rate will decrease ul¬ 
timately. The average collaboration rate of 75 actors and 
74 actresses show the mentioned pattern in the network, 
Fig. ([n|). Although the general behavior is similar for 
both curves, the mean value for number of collaboration 
for actresses reaches the maximum sooner. It behaves 
as if its dynamic equation has smaller order of differen¬ 
tiation, a, than the curve for actors collaboration. In 
other words, history effects and aging process have more 
severe impacts on actresses growth. Besides, studying 
the mean rate of collaboration on first and second half 
of the whole period shows no noticeable difference in dy¬ 
namical behavior. The above results are in remarkable 
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FIG. 6: Averaged number of collaboration per year for 75 
Oscar winner actors and 74 actresses. The rate of 
collaboration in both curves increases and after reaching a 
maximum value, it will decay gradually. 


agreement with findings of the numerical and simulations 
in the above sections as shown in Figure (l2|). 


VI. CONCLUSION 

Complexity and collective behavior are key character¬ 
istics of many natural and social systems. The past of a 
realistic systems can not be ignored in its dynamics and 
what is happening now. Even it could be said, since fu¬ 
ture expectations of a system are reselections of system 
goals, a complex system dynamic can be characterized 


by its goals too. Some mathematical models have been 
proposed to expand BA model by adding new terms to 
it’s differential equation, however, since these additional 
terms are not unique, there is no clear method to pre¬ 
fer one over the other. More over often these additional 
terms do not often have analytical solution. 

In attempt to describe how the history play role in the 
current events and growth dynamic of networks, we apply 
a kernel function as a average of the history on the stan¬ 
dard equation in BA model. This choice for the kernel 
shows interwoven of past events that leads to dynamic 
of the system. This approach leads to a fractional order 
BA differential equation governing network dynamic. 

According to results, this approach predicts a bound¬ 
ary for growth of members. Just by generalizing govern¬ 
ing equation in BA model to an equation with fractional 
order derivative, a more realistic dynamic for systems is 
achieved. As it can be seen in many real systems such as 
movie actors collaboration, networks of scientific papers, 
friendship networks and etc, degree of nodes progress for 
a certain period of time and then decay gradually after 
that. This change in topology of the system causes re¬ 
distribution in the members power. 

Moreover, even future goals of a system be it minimiz¬ 
ing system energy or adaptability with nature (in evo¬ 
lution) could shape its dynamic too. Future projection 
of kernel could be tough of as the way that future goal 
might effect one’s present decisions. Birds for example 
migrate seasonally to survive their life in next seasons. 
Therefore, considering the future goals could be an inter¬ 
esting study in network dynamics. 
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